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APPLICATION OF THE AMI PREDICTION METHOD 

TO A NUMBE|l OF AIRFOILS* 

F.A. Dvorak and B. Maskew 
Analytical Methods, Inc. 

SUMMARY 


A method for calculating the flow about airfoils up to and 
beyond the stall is described. It is an iterative procedure be~ 
tween potential flow and boundary layer solutions. The separated 
is modeled in the potential flow analysis using free vor- 
tex sheets which require an inner iteration to estcdjlish their 
shapes. The free vortex sheet length is an important parameter 
in the potential flow calculation. Results so far indicate a 
possible correlation between wake length and airfoil thickness/ 
^hord ratio. Calculated and experimental results are compared 
for a series of airfoils. 


INTRODUCTION 


Boundary layer separation is one of the least understood but 
most important of fluid flow phenomena affecting aerodynamic for- 
ces and moments. Its accurate modeling is essential to the esti- 
mation of airborne vehicle performance. Currently, reliance is 
placed on wind tunnel tests to determine the consequences of sepa- 
ration, a procedure which is not entirely free of doubt because 
of Reynolds number effects. Successful theoretical modeling of 
separation is limited to a small number of special cases, one of 
which is two-dimensional turbulent boundary layer separation from 
or diffusers. The first successful model for trailing- 
edge separation was developed by Jacob (ref. 1). With Jacob's 
model, the separation region is simulated using source fluid, the 
distribution of which is chosen to give constant pressure every- 
where in the separation region. In general, the method predicts 
the upstream pressure distribution in a satisfactory manner, al- 
though agreement with experiment for base pressure level is not 
consistent. 


M ..u ® separation model has been developed by Analytical 
Methods, Inc. which replaces the source distribution in the sepa- 
ration zone by a vortex wake model. This model is desc ribed in 

^ *Support was given by the U.S. Army Research Office, Research 
Triangle Park, N.C., for this work under Contract DAAG29-76-C-0019 . 
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some detail in reference 2, but is discussed herein for reasons 
of completeness. 

Symbols are defined in an appendix. 


SEPARATION MODEL 

Approximations for the Complete Flow Field 

An approximate model of the flow about an airfoi]. with a 
region of separation is shown in figure 1. It is assumed that: 

(i) The boundary layer anc free shear layers do not have signifi- 
cant thickness and, hence, can be represented as slip sur- 
faces; that is, streamlines across which there exists a jump 
in velocity. 

(ii) The wake region does not have significant vorticity and has 
constant total pressure (lower than the free-stream total 
pressure). It is, therefore, taken to be a potential flow 
region. 

The flow field in the potential flow is obtained using lin- 
early varying vortex singularities distributed on planar panels. 
The wake is represented by sheets of vorticity shed at the separa- 
tion points. 

The mathematical problem is to find the vorticity sheet 
strength such that the appropriate boundary conditions are met. 

The position of the vorticity sheet representing the free shear 
layer is not known a priori. 


Approximations For the Free Shear Layer 
(i) Wake Shape 

Initially, the streamlines are not known, and so the shapes 
of the free shear layers must be obtained iteratively starting 
from an initial assumption. Earlier calculations in which the 
vortex sheet shapes were obtained by iteration suggested an in- 
itial shape as follows. The upper and lower sheets are represen- 
ted by parabolic curves passing from the separation points to a 
common point downstream. The slope at the upstream end is the 
mean between the free stream direction and the local surface 
slope. (Indications from further calculations are that this 
starting slope should be streamwise for calculations beyond the 
stall.) Once the wake calculation begins, the initial slope and 
downstream position of each wake is determined by iteration. The 
final wake position represents the separating streamline. 


348 


I 


(II) Wake Length 

to that the results were sensitive 

to the lengtt of the free vortex sheets. Good correlation with 
experimental results was obtai-^ only with relatively short 

Such^a extending .ic to ,2c beyond the tral ling-edge. 

Such a model appears reasonable In the light of experimental ^ 
evidence! the separated wake does, in fait, close^iSJ^^down- 

as a result of the strong en train- 
la brought about by the rotation In the free shear 

with LpeJimen?^a"°® several comparisons 

function of the airfoil thickness tS cho?l Swo? 

This IS discussed in detail In reference 2. 

(III) Wake Pressure 

^jj® ^PP^o^i^ftlon of zero static pressure drop across the 
e shear layer is used to obtain an expression for the total 

Jons?d.r?!r S2e ?o?iex 

tSe SnoSr§y"““ velo- 
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since the vorticity, y„ - on the upper sheet. 

(The vorticity in the lower shear layer is v = v - v ^ 

'L inner outer 

The jump in total pressure across the shear layer is then 

~ ^’inner “ **outer 

l^inner 

- jPouter - ‘>p(v + Yu/2)‘ | 

= " pVYu - pvy^ 


given the boundary condition that the static pressure, p. 
jump in value across the shear layer. 


has no 
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Since the wake has constant total pressure (assumption (ii) ) / 
the jump in total pressure across the free shear layer is the 
same everywhere. 

Once the vorticity strengths of the individual panels rep- 
resenting the airfoil and of the vorticity sheets representing 
the wake are determined, the velocity at any point in the flow 
field can be calculated. 

The pressures are calculated from the velocities according 
tp the Bernoulli equation which is expressed non-dimensional ly 
as 



where C = - — ^ , q^= and AH = increase in total pres- 
sure over that”at infinity. Note that AH = 0 everywhere except 
in the wake region for which it was previously shown that AH 

= P^Yl* 


CALCULATION PROCEDURE 
Flow-Field Methods 


The overall calculation procedure is shown in figure 2, and 
involves two separate iteration loops. 


(i) Wake Shape Iteration .... , i.. 

The iteration loop for wake shape is the inner loop and in- 
volves the potential flow analysis only. Within this loop the 
separation points are fixed. The separation points may be loca- 
ted anywhere on a surface panel; they are not restricted to panel 

edge points. 

The wake shape is calculated as follows. Using the previous 
vorticity distribution, velocities are calculated at the panel 
mid-points on the free vortex sheets. The new "^ape is then 

determined by piecewise integration, starting at the _ separation 
points. The upper and lower sheet downstream end points, which 
were coincident in the initial wake, are allowed to 
ently in subsequent iterations. At each 

fluence coefficients at the surface control points are recalcula 
ted, and a new potential flow solution is obtained. 


The number of wake iterations is an input par^eter in the 
current version of the program; convergence criteria have not 
been investigated yet. 
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(ii) Viscous/Potential Flow Iteration 

Ai^ 4 - iteration loop takes the potential flow pressure 

over to the boundary layer analysis and returns with 
button^ points and with the boundary layer source distri- 

• e source distribution is determined directly from the 

boundary layer solution as a = A ^^ere is the stream- 

at the edge of the boundary layer, 
and 6 is the displacement thickness. The addition of this 
source distribution modifies the normal velocity, Vj^, at each 

ted^iegiSn^°^ point. The sources are set to zero in the separa- 
tion the new separa- 

tion points together with information from the previous iterated 

;nto potential flow solution is then obtained? a^dlfon. 
The outer iteration is terminated when the change in Co is below 
program^^"*^^ eight iterations is currently imposed within the 


Boundary Layer Methods 

boundary layer development on an arbitrarily- shaped two- 
dimensional lifting configuration with separated flow is very 
complex. A thorough and exact calculation of this development is 
domain of the time-dependent solution to the general 
Navicr Stokes equations. Unfortunately, the computer does not 
yet exist which is capable of handling such a problem in a 
reasonable time at a reasonable cost. Such a calculation is not, 
therefore, of practical interest to the aerodynamicist. Less 
ifficult or costly are the finite-difference boundary layer pro- 
grams now in existence. The amount of computer time required for 
each calculation still prohibits their use in an analysis proced- 
re of the type reported herein. Having made the above evalua- 
tion, one must conclude that if the objective is a viscosity- 
dependent calculation procedure of practical use to the aerodyna- 
micist for analysis, and, possibly for preliminary design, 

the method must be relatively simple to use and economic of com- 

achieved if integral boundary layer 
methods are used. In two dimensions, integral methods are typi- 
cally about 100 times faster than finite-difference methods. 

Tney can, however , be expected to break down in the region of 
separation where none of the boundary layer methods (including 

^® expected to be valid. It is anticipa- 
ued, therefore, that integral methods will suffice for most ap- 

plicatxons of interest to the aerodynamicist for Co predic- 
tion. 't'max 
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In 'those cases of special interest to the aerodynamicist^ 
such as the effect of area suction for boundary layer control or 
of roughness (rivets, etc.) on alternative boundary layer 

calculation modules are available. These methods are called as 
needed into the overall calculation procedure. A brief descrip- 
tion of the boundary layer methods is given in the following 
paragraphs . 

The laminar boundary layer development is calculated by 
Curie's method (ref. 4), an adaption of the well known method of 
Thwaites (ref. 5). The calculation proceeds ei'ther to laminar 
separation or to the end of the airfoil— —whichever occurs first. 
The calculated boundary layer development is then interrogated to 
determine if transition, laminar separation or forced transition 
(boundary layer tripping) has taken place. If any of these 
phenomena have occurred, the downstream flow is assumed to be 
turbulent. 

Methods for the calculation of turbulent boundary layers in 
two dimensions have been developed by many investigators . A 
review of these methods was made at a conference held in 1968 at 
Stanford University (ref. 6). One of the methods, an integral 
method by Nash and Hicks (ref. 7), compared very favorably with 
the more complex finite-difference methods. Now, several years 
later, the method remains an excellent approach for application 
to the current problem both in terms of accuracy and speed. 

If surface roughness or area suction are of interest, an 
alternate turbulent boundary layer method developed by Dvorak 
(refs. 8 and 9) can be called. This method is capable of predict* 
ing the downstream development and the skin friction drag of a 
turbulent boundary layer over a rough surface, or a surface with 
area suction boundary layer control. 

Turbulent boundary layer separation is predicted by either 
the Nash and Hicks or Dvorak methods when the calculated local 
skin friction coefficient reaches zero. 


DISCUSSION OF RESULTS 


The method was applied to a GA(W)-1 airfoil. This section 
shape represents a difficult test case and pressure distributions 
are available from experiments at NASA-Langley for a range of 
incidence. 

The first set of results, figures 3 through 5, are for a 
Reynolds number of 6.3 x 10® with a boundary layer trip at .08c. 
Figure 3 shows a very good agreement between the calculated and 
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experimental pressure distribution at 20 05 » mv 

culation took six viscous/potential flow it-i^T-a!?? *u 
three wake shaoe iteratinne flow iterations each with 

fiaurr® ”a!il ® incidence is shown in 

foil ItTRelnUrn^i^l 

9ood agreement with experiment up to c',^ , bS? Sle^^novSln 
the curve occurs 2 to 3 degrees later than in the experiment. 

. ShSin°on‘’tiJu?2s rfnd"8'^LrtL“MSi?rfof°thriif?^ 

(rS?°lo““irttrcasl of"S“LS“loog''airSiJ°“?l: 

Si S i°i: ‘k •*'® reparation reatta?hmen^c?it«iSi 

out the need lorTbl««”JS«“aSing‘Sf®^r''l'“"‘’ 

bubble bursting phenomenon!! ® laminar separation 


Comparisons between theory and experiment for rho Mf*. 

ar“fhoTO‘M^figSres*'l/lo*and^U?*' A °f Reynolds numbers 

and experimental variation with R^olL°num^r‘’u‘*shwn in 

ex^«Len;,a/c:rve'lrorre%«e^^^ 

sJS?ationS^bL''r® attempted, bit" •probiems^wttS'tilri^^L 
sults?^ bubble bursting criterion produced inconsistent re- 

calculations were made to demonstrate the cana- 

itiiiS!.°*Pi™rr?i''^h“ "r” * rangrif"igl2!o|!‘’* 

*“««">• The oorrespoidin?^pres- 

i«g ciii?L^‘i?L‘ i?e^srs(issfir2"Ji2 

■"•asurei'lif? ani1iar«e?n- 

ts given in the enclosed brackets. Figure 15 shows a 
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S''SACrSof2‘a!??oirl“fo%reer^^ dearUsiS^le-ol- 

attack. The agreement is surprisingly good. 


CONCLUSIONS 


v“?ie^; ofaSoUs over a range of Reynold, numbers,- 

(ii) post-stall behavior is best predicted for the trailing- 
edge type of stall; 


(iii) leading-edge and thin 
consx^ drably improved by a 
bursting criterion. 


airfoil stall prediction could be 
better laminar separation bubble 


(iv) the use of vortex sheets 
flow boundaries suggests that 
to unsteady flows. 


to represent the separated 
the model will be applicable 
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appendix 

Symbols 

airfoil chord 
drag coefficient 

lift coefficient 

maximum lift coefficient 

moment coefficient 

pressure coefficient 

total pressure 

static pressure 

free stream dynamic pressure 

Reynolds number, V c/v 

00 ' 

nonnsl v©locity 

average shear layer velocity 
free-stream velocity 

displacement thickness 

density 

vorticity 

boundary layer source strength = (V 6*) 
kinematic viscosity 
angle-of-attack 
upper surface 
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L lower surface 

oo value at free stream conditions 

e value at edge of boundary layer 


356 


f 


I 





REFERENCES 


1. Jacob, K. , "Computation of Separated, Incompressible Flow 
Around Airfoils and Determination of Maximum Lift", AVA Report 
66, A, 62, 1967. 

2. Maskew, B. and Dvorak, F.A.*, "The Prediction of C«„_„ Using a 
Separated Flow Model", Proceedings 33rd Annual Forum of the 
American Helicopter Society, Washington, D.C., May 1977. (To 
be published in the AHS Journal.) 

3. Seethar^, H.C. and Wentz, W.H., "Experimental Studies of Flow 
Separation and Stalling on a Two-Dimensional Airfoil at Low 
Speeds", NASA CR-2560 . 1975. 

4. Curie, H. , "A Two-Parameter Metliod for Calculating the Two- 

Dimensional Incompressible Laminar Boundary Layer", J.R. Aero. 
Soc. , voi. 71 ^ 1967. 

5. Thwaites, B., "Approximate Calculation of the Laminar Boundary 
Layer", Aero. Quart . , Vol. I, 1949. 

6. Kline, S.J., Morkovin, M.V., Sovran, G. and Cockrell, D.J., 
"Computation of Turbulent Boundary Layers", Proceedings 1968 
AFOSR-IFP Stanford Conference, Stanford University Press, 
Stanford, Calif., 1969. 

7. Nash, J.F. and Hicks, J.G., "An Integral Method Including the 
Effect of Upstream History on the Turbulent Shear Stress", 
Proceedings Computation of Turbulent Boundary Layers — 1968 
AFOSR-IFP Stanford Conference, Vol. 1, Stanford University 
Dept, of Mech. Eng., Stanford. Calif. 

8. Dvorak, F.A., "The Calculation of Turbulent Boundary Layers on 
Rough Surfaces in Pressure Gradient", AIaA J. , Vol. 7, No. 9. 
September 1969. 

9. Dvorak, F.A., "The Calculation of Compressible Turbulent Bound- 
ary Layers with Roughness and Heat Transfer", AIAA J., Vol. 10, 
No. 11, November 1972. 

10. McCullough, G.B. and Gault, D.E. , "Examples of Three Rept.;oen— 
tative Types of Airfoil Section ftall at Low Speed", NACA Tech. 
Note 2502, 1951. 

Pinkerton, R.N., "The Variation with Reynolds Number of Pres- 
sure Distribution over an Airfoil Section", NACA Report No. 

613, 1538. 


357 


REGION 1 - POTENTIAL FLOW REGION 

REGION 2 - BOUNDARY LAYER <hfpt<; 

REGION 3 - FREE SHEAR LAYER 

ppriOM 4 - WAKE REPRESENTING 

REGION 4 WAKE ^ FREE -SHEAR LAYERS 



REGION 3 


-VORTEX SHEET 
REPRESENTING 
BOUNDARY LAYER 

Flsure 1.- Mathematical flow model. 


£ 




Figure 2.- Calculation procedure. 
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Figure 5.- Comparisons of calculated and experimental lift and pitching- 
moment characteristics for the GA(W)-1 airfoil. Re * 6.3 ^ 10^. 
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Figure 6.“ Comparisons of calculated and experimental lift and pitching 
moment characteristics for the GA(W)-1 airfoil. Re = 2.1 10^. 




Figure 7.- Comparison of calculated and experimental lift characteristics 

for NACA 64A006 and 63009 airfoils. 
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Figure 8.- Comparison of calculated and experimental lift characteristics 
for NACA 63012 and 632018 airfoil sections. 
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Figure 9.- Comparison of calculated and experimental lift characteristics 
for an NACA 4412 airfoil. Re - 0.24 x iq 6 and 0.45 x 1 q6. 
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Figure 10.- Comparison of calculated and experimental lift characteristics 
for an NACA 4412 airfoil. Re - 0.9 x 106 and 1.8 x io6. 
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Figure 11.- Comparison of calculated and experimental lift characteristics 
for an NAioA 4412 airfoil. Re > 3.4 x 10^ and 6.3 x iq6. 
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Figure 13.- Calculated wake shape Tor an NACA 0012 airfoil at 90* incidence 
after six viscous-potential flow iterations, each with three wake-shape 
Iterations. 
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Figure lA.- Calculated pressure distribution on an NACA 0012 airfoil at 
90“ incidence. Re ■ 6.0 x 10^; Mach number, 0.2. 
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Figure 15.- Comparison of calculated and experimental lift characteristics 
for an NACA 0012 airfoil. Re * 6.0 x 10^; Mach number • 0.2. 
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